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Abstract. Spatially resolved, long-slit echellograms at 
different position angles of the bright, low excitation plan- 
etary nebula NGC 40 indicate that, the higher is the gas 
excitation, the faster is the radial motion, thus confirm- 
ing the overturn of the Wilson law already suggested by 
Sabbadin and Hamzaoglu (1982). 

New reduction procedures, giving the radial trends of 
the electron density and of the ionic and chemical abun- 
dances, were applied to NGC 40; they show that: 

- the radial matter distribution has a sharp "bell" pro- 
file with peaks up to 4000 cm" 3 ; 

- the ionization structure is peculiar, indicating the 
presence of chemical composition gradients within the 
nebula: the innermost regions, hydrogen depleted, are es- 
sentially constituted of photospheric material ejected at 
high velocity by the WC8 nucleus. 

Moreover, detailed H + , ++ and N + tomographic 
maps, giving the spatial ionic distributions at four position 
angles, are presented and discussed within the interacting 
winds evolutionary model. 
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1. Introduction 

NGC 40 (PNG120.0+09.8, Acker et al., 1992; Figure 1) is 
a very low excitation planetary nebula powered by a WC8 
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Fig. 1. Broad-band R image of NGC 40 (exposure time 
180s, seeing 0.85") obtained with the 3.5m Italian Na- 
tional Telescope (TNG), revealing the intricate Ha+[NII] 
nebular structure. North is up and East is on the left. 



star presenting a mass-loss rate of the order of 10 6 - 10 8 
M yr _1 with wind velocities between 1800 and 2370 km 
s _1 (Cerruti-Sola and Perinotto, 1985; Bianchi, 1992). 

The nebula is described by Curtis (1918) as a trun- 
cated ring, from the end of which extend much fainter 
wisps; the brighter central portion is 38x35" in PA=14°, 
while the total length along this axis is about 60". Deep, 
narrow-band imagery by Balick et al. (1992) reveals the 
presence of an external network of knotty floccules and 
smooth filaments. Meaburn et al. (1996) identify two 
haloes around the barrel-shaped core, and a jet-like fea- 
ture projecting from it. 

Detailed low resolution spectroscopy (Aller and 
Czyzak, 1979, 1983; Clegg et al., 1983) indicates point to 
point changes of the excitation and chemical abundances 
which are typical of planetary nebulae, without apparent 
contamination by the fast stellar wind. 

A kinematical analysis secured by Sabbadin and 
Hamzaoglu (1982) suggests the presence of peculiar mo- 
tions in the central part of NGC 40, with the expansion 
velocity larger in [OIII] than in Ha and [Nil]. 

More recently, Meaburn et al. (1996) obtained long-slit 
Ha and [Nil] high resolution spectra of the bright core, of 
the two haloes and of the jet-like feature. These authors 
showed that the outer filamentary halo is practically in- 
ert and that the motion of the inner, diffuse, spherical 
halo mimics the one observed over the bright nebular core; 
moreover, the jet-like structure is kinematically associated 
with the receding end of the barrel shaped-core and does 
not present any of the characteristics expected of a true 
jet. 

In order to study in detail the nebular physical condi- 
tions and to apply an original procedure giving the spatial 
matter distribution along the cross-section covered by the 
spectroscopic slit, on December 1998 we obtained a num- 
ber of echellograms at different position angles of a dozen 
of winter planetary nebulae (NGC 40, 650-1, 1501, 1535, 
2022, 2371-2, 2392, 7354 and 7662, J 320 and 900, A 12 
and M 1-7). 

In this paper we report the results derived for NGC 
40. In Section 2 we present the observational material; in 
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Section 3 we discuss the expansion velocity field; in Sec- 
tion 4, the electron temperature and the electron density 
are analysed and a new method to determine the radial 
distribution of the electron density in expanding nebulae 
is presented. The radial ionization structure of NGC 40 is 
given in Section 5 and the resulting chemical composition 
gradients in Section 6. Section 7 shows the H + , ++ and 
N + tomographic maps derived at four position angles; a 
short discussion is contained in Section 8 and conclusions 
are drawn in Section 9. 

2. Observational material 

Broad-band U, V and R imagery of NGC 40 was obtained 
with the Optical Imager Galileo (OIG) camera mounted 
on the Nasmyth A adapter interface (scale=5.36 " mm -1 ) 
of the 3.5m Italian National Telescope (TNG, Roque de los 
Muchachos, La Palma, Canary Islands). The OIG camera 
was equipped with a mosaic of two thinned and back- 
illuminated EEV42-80 CCDs with 2048x4096 pixels each 
(pixel size=13.5 microns; pixel scale= 0.072 " pix -1 ). 

Though these images were taken during the testing 
period of the instrument (when interference filters were 
unavailable) , they reveal new morphological details of the 
nebula supporting our spectroscopic results and will be 
presented later-on (Section 5). 

Spatially resolved, long-slit spectra of the bright core 
of NGC 40 in the range AA4500-8000 A (+flat fields+Th- 
Ar calibrations+ comparison star spectra) were obtained 
with the Echelle spectrograph (Sabbadin and Hamzaoglu, 
1981, 1982) attached to the Cassegrain focus of the 182cm 
Asiago telescope, combined with a Thompson 1024x1024 
pixels CCD. 

Four position angles were selected: 20° (apparent ma- 
jor axis), 110° (apparent minor axis), 65° and 155° (in- 
termediate positions). In all cases the slit-width was 0.200 
mm (2.5" on the sky), corresponding to a spectral resolu- 
tion of 13.5 km s _1 (1.5 pixel). During the exposures the 
slit grazed the bright central star (we tried to avoid most 
of the continuum contamination; the faint stellar spectrum 
present in our echellograms was used as reference for the 
nebular centre and to correct the observations for seeing 
and guiding errors). 

All spectra were calibrated in wavelength and flux in 
the standard way using the IRAF data reduction software. 

The following nebular emissions were detected: H/3, 
A4959 A and A5007 A of [OUT], A5755 A of [Nil], A5876 
A of Hel, A6300 A and A6363 A of [OI] (the first at the 
order edge; both partially blended with night-sky lines), 
A6548 A and A6584 A of [Nil], Ha, A6578 A of CII, A6717 
A and A6731 A of [SIT], A7135 A of [Arlll], A7320 A and 
A7330 A of [Oil] (at the extreme edge of the order; out of 
focus) and A3726 A and A3729 A of [Oil] (second order). 

All line intensities (but not the second order doublet 
of [Oil]) were de-reddened by fitting the observed Ha/H/3 
ratio to the one computed by Brocklehurst (1971) for 



Te=10 4 K and Ne=10 4 cm 3 . We derive a logarithmic ex- 
tinction at H/3, c(H/3)=0.60 ± 0.10, to be compared with 
the values of 0.28, 0.65 and 0.70 obtained by Kaler (1976), 
Aller and Czyzak (1979) and Clegg et al. (1983), respec- 
tively. The fair agreement (within 10%) between our inte- 
grated line intensities and those reported in the literature 
from low resolution spectra (see, for example, Clegg et al., 
1983) induced us to adopt I(3726+3729)/I(Ha)= 1.48, as 
reported by these authors. 

A large variety of emission structures is present in 
NGC 40; for illustrative purposes, some examples are 
given in Figure 2. In these reproductions the faint stel- 
lar continuum was removed and the observed intensities 
enhanced (by the factor given in parenthesis) to make each 
line comparable with Ha. 

Our spectra at PA=110° (apparent minor axis, Fig. 
2a) were centred 2" South of the central star. All emissions 
present the classical bowed shape expected for expanding 
shells of different mean radii r and thicknesses t, ranging 
from r = 16" for [OIII] and CII to r = 19" for [Nil] and 
[SII] and from t = 0.20r for [Nil] and [SII] to t = 0.25r for 
[OIII] and CII (as also suggested by the integral intensity 
distributions along the slit). 

At PA=20° (apparent major axis; slit centre 2" West 
of the star; Fig. 2b) low excitation lines (such as the [SII] 
red doublet) show knotty structures and fail to close at 
either polar cap. High excitation emissions (for instance 
[OIII]) appear more homogeneous; they merge in a diffuse 
zone internal to the polar caps; the smooth intensity dis- 
tribution and the large velocity spread suggest that [OIII] 
completely fills this part of nebula; moreover, the faintness 
(or absence) of an external low excitation counterpart in- 
dicates that in these directions the main body of NGC 40 
is density bounded rather than ionization bounded. 

Spectra taken al PA=65° (Fig. 2c) were centred 2" S- 
E of the central star; they appear quite regular, but for 
a deformation of the blue-shifted component in the E-NE 
quadrant. 

At PA=155° (centre of the slit 2" S-W of the star; Fig. 
2d) low excitation lines (such as the [Nil] line at A6584 A) 
merge in the N-NW sector, in correspondence of a bright 
knot. The S-SE region is strongly depauperated of gas and 
only high excitation emissions (like the extremely faint 
CII line at A6578 A) present a close structure; as already 
seen for PA=20°, the weakness of external low excitation 
emissions suggests that in these directions the nebula is 
density bounded. 

3. The expansion velocity field 

First of all we measured the peak separation at the centre 
of each nebular emission. Results are presented in Table 
1, where ions are put in order of increasing ionization po- 
tential. 

We find a complete overturn of the Wilson law: in NGC 
40 the higher is the ionization, the faster is the motion 
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Fig. 2. Representative examples of emission line structures observed in NGC 40 at different position angles. Each line 
intensity was multiplied by the factor given in parenthesis to make it comparable with Ha. 
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Fig. 3. Observed (left panel) and deprojectcd (right 
panel) expansion velocity field at PA=110° (apparent mi- 
nor axis). Triangles are for [Nil] and circles for [OIII]. 

Table 1. Expansion velocities 



A 

(A) 


ion 


LP. 

(eV) 


2 Vea;p 

(km s" 1 ) 


notes 


6300 


[OI] 





51 ±2 


1 


6717-6731 


[SII] 


10.4 


51 ±2 




6563 


HI 


13.6 


50 ±1 




3726-3729 


[Oil] 


13.6 


54 ±2 


2 


6584 


[Nil] 


14.5 


51 ±1 




6578 


CII 


24.4 


53 ±2 




5876 


Hel 


24.6 


53 ±2 




7135 


[Arlll] 


27.6 


54 ±2 




5007 


[OIII] 


35.1 


57 ±1 





1 = partially blended with night-sky emission 

2 = second order lines 

(thus confirming the suggestion by Sabbadin and Hamza- 
oglu, 1982, of a peculiar [OIII] expansion velocity). 

A further peculiarity is represented by [Oil] , which ex- 
pands faster than the other low ionization species. This 
[Oil] velocity-excess is related to the very low degree of 



Fig. 4. Same as Fig. 3, but for PA=20° (apparent ma- 
jor axis). Glossary: N=main nebula, P=polar caps (blue- 
shifted at south, red-shifted at north), M=extremely faint 
Ha and [Nil] mustaches (red-shifted at south, blue-shifted 
at north) emerging from the main nebular body. 

excitation of the nebula: since the + zone stops at 35.1 
eV (while, for example, S + falls at 23.4 eV and N+ at 
29.6 eV), it extends well inwards, as also indicated by the 
[OIII] weakness. So, in NGC 40, [Oil] must be considered 
a medium-high excitation species (at least for the dynam- 
ical point of view). 

Clearly, the expansion velocities reported in Table 1 
are only mean values and do not describe the complex 
kinematics of the nebula. To better illustrate this, we have 
selected two position angles (110°=apparent minor axis; 
20°=apparent major axis) and chosen [Nil] for low excita- 
tion regions and [OIII] for the high excitation ones (" high" 
at least for NGC 40!). 

The regular ellipses shown by these ions at PA=110° 
(Figure 3, left panel) are typical of expanding shells; the 
absence of tilt indicates that the nebular cross-sections 
are either circular, or elliptical (in this second case we are 
aligned with one of the axes). In fact, the de-projection 
of the observed expansion velocities (assuming a simple, 
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Fig. 5. [SII] electron density trends observed at PA=110° 
for the blue-shifted (dotted curve) and red-shifted (dashed 
curve) peaks and for the integrated spectrum (solid curve). 
"Position a" of Clegg et al. (1983) is indicated. 

direct position-speed correlation) gives the "true" expan- 
sion velocities presented in the right panel of Figure 3, i.e. 
a constant value for each ion. Similar trends are obtained 
also at PA=65° and in the N-NW sector of PA=155° . 

The [Nil] and [OUT] expansion velocity fields observed 
at PA=20° (major axis) are given in Figure 4 (left panel). 
Here the situation is more complicated, mainly with re- 
spect to [OIII]. In fact, while [Nil] presents sharp emis- 
sions along the slit, [OIII] describes an ellipse which is 
well defined in the main part of the nebula, but be- 
comes smooth and broad to the north and south. The de- 
projection of the observed velocities (right panel of Fig- 
ure 4) points out the acceleration suffered by the polar 
(less dense) nebular material. A similar behavior occurs 
at PA=155°, in the S-SE sector. 

For a better understanding of the peculiar kinematics 
of NGC 40, a detailed analysis of the nebular physical 
conditions can be useful. 

4. Electron temperature (Te) and electron 
density (Ne) 

The only Te diagnostic present in our spectra is the 
line intensity ratio 6584/5755 A of [Nil]. Unfortu- 
nately, the auroral emission is too weak for a detailed 
analysis; for the brightest parts of NGC 40 we ob- 
tain I(6584)/I(5755)=125±20, corresponding to Te=7900 
±200 K (assuming Ne=10 4 cm -3 ). Thus, according to 
Aller and Czyzak (1979, 1983) and Clegg et al. (1983), 



Fig. 6. Same as Fig. 5, but for PA=155°. "Position b" of 
Clegg et al. (1983) is indicated. 

we will adopt Te=8000 K for low ionization regions (O , 
N+, S+, 0+) and Te=10000 K for the others (0++, Hc+, 
C+ Ar++). 

Our iVe study is essentially based on [SII] lines at 
A6717 A and A6731 A. 

To test the reliability of our measurements we can 
make a comparison with the results of Clegg et al. (1983), 
who observed at low spectral resolution a patch of neb- 
ulosity located at north-west of the central star. At "po- 
sition a" they found I(6717)/I(6731)=0.68, corresponding 
to Ae=2200 cm" 3 (for Te=8000 K). "Position a" practi- 
cally coincides with the bright western edge of our spectra 
taken at PA=110° (apparent minor axis). 

The [SII] density distribution we observed at this posi- 
tion angle is shown in Figure 5 for the blue-shifted and the 
red-shifted peaks, and also for the integrated, low resolu- 
tion spectrum (we have averaged our echellograms over 
a velocity range of 312 km s , which corresponds to 
the spectral resolution used by Clegg et al., 1983). The 
estimated Ne accuracy varies from ±15% for high den- 
sity peaks (generally -but not always- coinciding with the 
strongest emissions), to ±50% for deep valleys (weakest 
components). 

For "position a" the integrated spectrum gives 
Ne =2100 cm~ 3 , in good agreement with Clegg et al. 
(1983), but high resolution measurements indicate a den- 
sity of 3100 cm~ 3 . Note that in Figure 5 the maximum 
density of the integrated spectrum (3000 cm -3 ) is reached 
in a relatively weak knot projected onto the central star; 
it corresponds to a peak (3200 cm -3 ) in the blue-shifted 
component. 
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Fig. 7. Relative (left ordinate scale) and absolute (right 
ordinate scale) radial electron density distributions at 
PA=110° from the Ha (solid curve) and [Nil] (dashed 
curve) zero-velocity columns. 

As a further example, in Figure 6 we present the [SII] 
density distribution observed at PA=155°; symbols are 
as in Figure 5. In this case our N-NW peak is partially 
superimposed to "position b" of Clegg et al. (1983). Once 
again their electron density (A^e=2100 cm~ 3 ) coincides 
with our integrated value, but is considerably lower than 
the peak measured at high resolution (A^e=2900 cm -3 ). 

CAVE AT '-.strictly speaking, our Ne estimates too 
must be considered as lower limits, due to ionization ef- 
fects. In fact, as we will see in the next section, the com- 
plete recombination of S ++ occurs externally to the top of 
the density distribution; this implies that the peak of the 
[SII] intensity distribution is slightly outwards with respect 
to the density maximum, leading to a small underestimate 
of Ne derived from the [SII] line intensity ratio. 

Besides the density peak along the cross-section of the 
nebula, our high resolution spectra allow to derive the 
radial distribution of the electron density using the " zero- 
velocity column" of each emission line; this is the central, 
rest-frame pixel column containing the ionized gas moving 
in the neighborhood of -30 km s _1 (mean radial velocity 
of the whole nebula). In practice, the zero- velocity column 
contains the material which is expanding normally to the 
line of sight; it is important because it insulates a defi- 
nite slide of nebula which is unaffected by the expansion 
velocity field. 

The intensity distribution observed in the zero- velocity 
column of each nebular emission must be cleaned of the 
broadening effects due to: 




u i i i i i i i i i i i i 

40 20 20 

arcsec 



Fig. 8. Radial electron density profiles obtained from the 
Ha zero-velocity columns at PA=20°, 65° and 155°. 

- instrumental profile, corresponding to a Gaussian 
having FWHM=13.5 km s _1 ; 

- thermal motions; for Te=8000 K they amount to 19 
km s _1 for H, 9.5 km s _1 for He and 5 km s _1 (or less) 
for the heavier elements; 

- turbulent motions, which affect all lines in the same 
way; following Meaburn et al. (1996) in NGC 40 they are 
6-8 km s _1 ; 

- fine structure of Ha and A5876 A of Hcl. Following 
Dyson and Meaburn (1971) and Dopita (1972), the seven 
Ha components can be shaped as the sum of two equal 
Gaussian profiles separated by 0.14 A. A5876 A of Hel 
has six components: five of them are concentrated in the 
narrow range AA5875.5987 - 5875.6403 A (they were con- 
sidered a single emission at A5875.62 A); the sixth one, 
located at A5875.9663 A, is 4.4 times weaker (and then it 
was neglected). 

We have derived the corrected intensity of the zero- 
velocity columns by deconvolving the observed profile at 
any given nebular position along the slit into a series of 
Gaussians spaced of one pixel along the dispersion and 
having: FWHM=25.8 km s^ 1 for Ha, 17.9 km s" 1 for 
A5876 A of Hel and FWHM=16.0 km s" 1 for the other 
lines. We have cautiously stopped the computation when 
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Iobs=2I CO rr, corresponding to an uncertainty in Ne of 
±30%. 

The last step is to deconvolve the corrected zero- 
velocity columns for guiding errors and seeing by means of 
the profile, along the slit, of the stellar continuum present 
in the spectrum. 

At each position of the zero-velocity column the in- 
tensity is proportional to NeNi and, in the case of total 
ionization and Te=constant, IocAe 2 . This is valid, sepa- 
rately, for all ions. 

As an example, the Ha and A6584 A of [Nil] zero- 
velocity columns derived at PA=110° (apparent minor 
axis) are presented in Figure 7; densities are normalized 
to the strongest peak (left ordinate scale). As expected, 
A^e(H + ) and Ne(N + ) practically coincide in the external 
parts (up to the maxima) , and tend to diverge moving in- 
wards, due to the gradual double ionization of nitrogen. 
Let's consider the H+ curve; if we scale the western peak 
to Ne—3100 cm~ 3 , as previously derived from the [SII] 
line intensity ratio (see Fig. 5), we obtain the true density 
distribution along the central, zero- velocity column (right 
ordinate scale). To be noticed that the density of the east- 
ern peak results to be 2600 cm~ 3 , in fairly good agreement 
with the value derived from [SII] diagnostics. 

In the same way we obtained the radial density trends 
at the other position angles; they are given in Figure 8. 
Note that the profiles observed at PA=20° correspond 
to the northern (Ae=1000 cm~ 3 ) and southern (Ae=600 
cm~ 3 ) polar caps, and that only a weak emission (Ae=400 
cm~ 3 ) is visible in the S-SE sector of PA=155°. 

In the end, the "bell" profile seems to be a general char- 
acteristic of the radial matter distribution within NGC 40. 

A final note on the Ne problem in this nebula concerns 
[Oil] . We have repeated the previous procedures using the 
3726/3729 diagnostic ratio; the [Oil] density trends are, 
in all cases, very close to the [SII] ones, but the actual 
density is systematically lower by a factor of almost two. 
To explain this curious feature, we must examine in detail 
the ionization structure of NGC 40. 

5. Radial ionization structure 

The detailed radial ionization structure of an expand- 
ing nebula can be derived by comparing the zero- velocity 
columns of different lines. 

To obtain the ionic abundances (relative to H + ) of 0°, 
0+, 0++, N+, C+, S+, He+ and Ar++ we have applied 
the classical procedure (see, for instance, Peimbert and 
Torres-Peimbert, 1971; Barker, 1978; Aller and Czyzak, 
1979, 1983). 

For all ions, but not [Arlll], we have used the same 
atomic constants adopted by Clegg et al. (1983) in their 
low resolution study of NGC 40; for [Arlll], not mentioned 
by these authors, transition probabilities were taken from 
Mendoza (1983) and collision strengths from Galavis et 
al. (1995). 
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Fig. 9. Radial ionization structure of NGC 40 at PA=155° 
obtained from the zero- velocity columns. Details are given 
in the text. 



Since the H+ central column is practically absent in 
the S-SE sector of PA=155° (while other ions are present) , 
only lower limits to the relative ionic abundances can be 
obtained; the same situation occurs at PA=20° (major 
axis), internally to the polar caps. 

In all the remaining cases (where the radial ioniza- 
tion structure of the nebula can be studied in detail) most 
ions present the expected trends: medium-high ionization 
species (0 ++ , Ar ++ , He+) rapidly decrease outwards; the 
opposite for low ionization species, such as O , N + and 
S + . Once again, + shows an anomalous behavior (also 
C + has similar characteristics, but the CII A6578 A line 
is so weak that a quantitative analysis is hazardous) . 

The radial ionization structure of NGC 40 is synthe- 
sized in Figure 9 for PA=155°; this position angle is a 
perfect compendium of the different situations occurring 
in the nebula, since it presents a well defined H + structure 
in the N-NW part and practically no hydrogen emission 
at S-SE (see Figure 8). Moreover, its N-NW peak corre- 
sponds to "position b" of Clegg et al. (1983), so that a 
direct comparison with their results can be made. 

Let's first consider the N-NW sector of Figure 9; here 
both + /H + and ++ /H + increase inwards (while we 
would expect opposite trends, since the + ^0 ++ equi- 
librium shifts on the left going outwards, so that when 
[Oil] increases [OIII] falls, and vice versa). Note also that 
the ++ curve lies well below the + one, indicating that 
the stellar flux is unable to completely ionize + . Ioniza- 
tion structures similar (qualitatively and quantitatively) 
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to the one observed in the N-NW sector of PA=155° were 
obtained at PA=65° and 110°. 

Look now at the left part of Figure 9. In the S-SE 
sector of PA=155° the H + zero-velocity column shows 
only a weak emission at about 30" from the center (and 
here the N+/H+ trend was derived), whereas high ion- 
ization species are present also at intermediate positions, 
giving the lower limits to the He+/H+, C+/H+, 0+/H+ 
and ++ /H + abundances indicated in the figure. The de- 
gree of excitation in these hydrogen-depleted zones re- 
mains low, being N(0 + )/N(0 ++ )^5. Similar low density, 
hydrogen-depleted regions are present at PA=20° (major 
axis), internally to the polar caps. 

An independent confirmation of the existence of a pe- 
culiar ionization structure in NGC 40 comes from the 
broad-band U, V and R imagery obtained with the 3.5m 
Italian National Telescope (TNG). 

The R band (Figure 10 upper left) is dominated by Ha 
and [Nil], with a small contamination of [SII] at A6717 A 
and A6731 A. 

Over 90% of the U light (Figure 10 upper right) comes 
from the [Oil] doublet at A3726 A and A3729 A; the re- 
maining is due to high Balmer lines. 

Finally, in the V band (Figure 10 lower left) H/3 con- 
tributes for 55% and the [OIII] doublet at A4959 A and 
A5007 A for 45%. 

Since the H/3 morphology of this low excitation plan- 
etary nebula is very similar to the Ha+[NII] one, we can 
extract the [OIII] component of the V image by subtract- 
ing the R image (appropriately scaled in intensity) to the 
V image. The resulting [OIII] appearance of NGC 40, 
shown in Figure 10 lower right, noticeably differs from 
the Ha+[NII] one; this morphological peculiarity, known 
since Louise (1981), stands out in great detail (greater 
than before; see, for instance, Balick, 1987 and Meaburn 
et al., 1996). 

The apparent [Oil] /[OIII] distribution over the neb- 
ula can be obtained from the ratio U/[OIII] and 
the(Ha+[NII])/[OII] distribution from R/U; results are 
shown in Figures 11. 

The large stratification effects present in these 
last maps indicate that all the three intensity ratios 
[Oil] /[OIII], (Ha+[NII])/[OII] and (Ha+ [Nil])/ [OIII] de- 
crease inwards, i.e. all the three ionic ratios + /0 ++ , 
(H++N+)/0+ and (H++N+)/0++ decrease inwards. In 
conclusion: all the three ionic ratios ++ /0 + , + /H + 
and ++ /H + increase inwards (being the H + distribution 
very similar to the N+ one). This is in perfect qualitative 
agreement with the results obtained from high resolution 
spectrography. 

The obvious conclusions emerging from the peculiar- 
ities in the ionization structure of NGC 40 are drawn in 
the following, short section. 



Fig. 10. Broad-band imagery of NGC 40 (seeing 0.85") 
obtained with the 3.5m Italian National Telescope (TNG). 
The R image (upper left, exposure time 180s) is dominated 
by Ha and [Nil] and the U image (upper right, 600s) by 
the [Oil] doublet at A3726 and A3729 A. The V image 
(lower left, 180s) is the sum of H/3 (55%) and [OIII] at 
A4959 and A5007 A (45%). The lower right panel shows 
the [OIII] appearance of NGC 40; it was obtained by sub- 
tracting the R image (appropriately scaled in intensity) to 
the V image. 



Fig. 11. Left panel: apparent [Oil] /[OIII] distribution in 
NGC 40. [OIII] is weak in the outermost strata (white), 
it gradually increases inwards, reaching its maximum 
(black) internally to the polar caps. Right panel: apparent 
(Ha+[NII])/[OII] distribution in NGC 40. This intensity 
ratio presents a maximum (white) in the external, bright 
shell, it decreases inwards, and has a minimum (black) in- 
ternally to the polar caps. Same orientation and scale as 
Fig. 10. 

6. Abundance gradients towards NGC 40. 

We believe that the radial ionization structure of NGC 
40 indicates the presence of abundance gradients across 
the nebula (although we cannot exclude a minor effect 
due to electron temperature variations). Take oxygen, for 
instance. 

The general expression 0=0°+0 + +0 ++ +ctc. in the 
case of NGC 40 reduces to 0=0 + , i.e. the observed 
+ /H + ratio is a good match of the total oxygen abun- 
dance. We obtain 0/H=6±l xl0~ 4 in the main, denser 
regions (to be compared with the values of 6.0xl0 -4 and 
8.4xl0 -4 given by Aller and Czyzak, 1983, and Clegg et 
al., 1983, respectively). O/H rapidly increases inwards; at 
the observational limit it is larger than 5xl0~ 3 . 

Less accurate abundance gradients can be derived in 
a similar manner for C/H (lxl0~ 3 to >3xl0~ 2 ) and for 
He/H (4xl0~ 2 to >lxl0- 1 ). 

It is very tempting to attribute these abundance trends 
to contamination of the stellar wind (we remember that 
detailed analyses of the optical spectrum of the WC8 nu- 
cleus carried out by Leuenhagen et al., 1996, give an upper 
limit for the mass fraction of hydrogen of 2%: moreover, 
Hc/C=0.8 and C/O— 5). In this scenario, our internal, low 
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density, hydrogen-depicted regions are essentially consti- 
tuted of enriched material recently emitted by the central 
star. Moving outwards, this gas gradually mixes with the 
shell's gas and abundances rapidly fall to the nebular val- 
ues. 

The general expression locNeNi (valid in the zero ve- 
locity column) allows to estimate the electron density in 
the innermost, fast moving, hydrogen-depleted regions. An 
indicative value of 1 - 2 cm" 3 is obtained, corresponding 
to a few 10~ 24 g cm~ 3 , and (for a distance of 1100 pc) to 
a matter flux of a few 10 -8 M© yr _1 . This must be the 
order of magnitude of the present central star's mass-loss. 

7. Tomography 

In order to obtain the spatial matter distribution within 
a planetary nebula we will analyse the observed emission 
line structures by means of an iterative procedure which 
is the implementation of the method originally proposed 
by Sabbadin et al. (1985, 1987 and references therein) for 
plate echellograms. In this approach, the relative density 
distribution of an emitting region along the cross-section 
of nebula covered by the slit is obtainable from the radial 
velocity, FWHM and intensity profile. 

The calibration to absolute densities is straightforward 
if [SII] or [Oil] doublets are present (as in the case of NGC 
40); for density bounded planetary nebulae other diagnos- 
tics (such as 5517/5537 of [C1III] or 4711/4740 of [ArlV]) 
can be used but, due to the intrinsic weakness of these 
lines, results are often uncertain; so, for high excitation 
planetary nebulae, we prefer to adopt the surface bright- 
ness method (it will be described in a forthcoming paper 
dedicated to NGC 1501). 

The tomographic analysis contains a large amount of 
physical informations; for reasons of space too, this first 
application to NGC 40 intends to be a concise test of re- 
liability (i.e. a semi-quantitative comparison with the ob- 
servational results independently accumulated in the pre- 
vious sections). 

Amongst the different maps (intensity, density, frac- 
tional ionization, ionic abundance etc.) obtainable with 
our reconstruction method, we have selected ^/NeNi; H+ 
was chosen as ionization marker, and N + and ++ to 
symbolize low and high excitation zones, respectively. 

Since for H + we have that ^/NeNi is proportional to 
Ne, the H + maps give the relative density distributions 
along the cross-sections of NGC 40 covered by the slit; 
this is roughly true also for N+ (its ionization occurring 
only in the innermost parts), but not for ++ , due to the 
incomplete + ionization. 

The gray-scale ^/NeNi maps derived for H+, ++ and 
N+ at PA=110°, 20°, 65° and 155° are shown in Figure 
12. 

At PA=110° (apparent minor axis) H + describes a cir- 
cular, knotty ring; N + closely mimics H + (but, obviously, 
it is a bit sharper), while ++ is emitted in the internal 



regions and rapidly drops when Ne raises. To be noticed 
the opposite behaviors of H + and ++ , in the sense that 
bright H+ zones tend to be associated with faint ++ 
regions, and vice versa; the patch studied by Clegg et al. 
(1983, their "position a") corresponds to our western edge. 

At PA=20° (apparent major axis) the H+ and N+ ring- 
structures appear disrupted at north and south, where the 
two polar caps are formed; unfortunately, in these repro- 
ductions we have lost the two extremely faint (and fast 
moving) low excitation mustaches attached to the main 
nebular body (blue-shifted at north, red-shifted at south). 
The elliptical structure shown at this position angle by 
++ testifies the existence of hydrogen-depleted zones in- 
ternally to the polar caps (which, in our opinion, are es- 
sentially constituted of stellar wind). 

The Ne map derived from H + at PA=65° describes a 
distorted ring presenting a deformation in the approaching 
gas of the E-NE quadrant; the bright knots visible in both 
the H + and N + maps completely disappear in the ++ 
one, which presents a more homogeneous distribution. 

Finally, at PA=155° the ring structure shown by H+ 
and N + in the N-NW sector is completely absent at S- 
SE; once again, ++ indicates the presence of hydrogen- 
depleted gas within the S-SE region; the N-NW border 
corresponds to "position b" of Clegg et al. (1983). 

Quantitative maps, giving the true electron density 
distribution at the four position angles, can be obtained 
by calibrating the H + maps of Fig. 12 with the [SII] line 
intensity ratio. Results are shown in Figure 13, where the 
nebular centre corresponds to the position of the exciting 
star (which, actually, lies 2" over or below the plane of the 
figure). 

8. Discussion 

The general rule given by Olin Wilson (1950) in his clas- 
sical work on the expansion velocity of planetary nebulae 
-"high-excitation particles show smaller separation, and 
low-excitation particles higher separation"- is definitely 
violated in NGC 40, which presents a positive inwards ve- 
locity gradient. Moreover, this velocity gradient is small 
enough (AV ~ 3 — 4 km s _1 for the main shell) to rouse 
the suspicion that some nebular lines can be produced by 
resonance scattering. 

The problem has been analysed in detail by Clegg et 
al. (1983), who tried to explain the abnormal CIV A1549 
A intensity observed in NGC 40 in terms of stellar light 
scattered by the nebula. Assuming a spread in velocity of 
20 km s _1 , these authors obtained that the stellar flux 
at A1549 A is nine times larger than the A1549 A nebular 
scattered flux and concluded that scattering is inadequate, 
but pointed out that "this mechanism will be most effec- 
tive if velocity gradients are small" . Decidedly, our new 
AV value re-opens the question. 

We have seen that the peculiarities found in the ioniza- 
tion structure strongly suggest the presence of abundance 
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Fig. 12. Gray-scale ^jNeNi sp atial m aps of NGC 40 in H + (left), ++ (center) and N + (right) at different position 
angles. For H + we have that \J NeNi is proportional to Ne; this is approximately valid also for N + , due to the very 
low degree of excitation of the nebula, but not for ++ , since the + ionization is incomplete . 



Fig. 13. Isophotal density maps of NGC 40 obtained from H + . Contour levels: dashed curve=250 cm 3 , solid 
curves=500, 1000, 1500, 2000, 2500, 3000 and 3500 cm~ 3 . Same orientation as Fig. 12. In each map the central 
star is located 2" over or below the plane of the figure. 



gradients within NGC 40; the responsible of these gradi- 
ents is the fast stellar wind. This is not a surprise since, 
following Bianchi and Grewing (1987), the mass of the 
central star is close to the C-0 core of the progenitor star. 

It is interesting to note the discrepancy found by these 
authors between the stellar temperature determined from 
IUE observations (90000 K) and the much lower value of 
30000 K inferred from the nebular ionization. They over- 
came the impasse with a carbon curtain screening the high 
energy photons; this 2 - 5xl0 18 cm~ 2 CII column density 
should be located at the inner edge of the shell, where the 
fast wind from the nucleus interacts with the nebular gas. 

The hypothesis is suggestive, but some doubts remain, 
essentially for lack of space. In fact, on one side our ion- 
ization maps indicate that the CII layer producing the 
requested column density must have a thickness of the 
order of 10 19 cm, on the other side the mean nebular ra- 
dius of NGC 40 (assuming a distance of 1100 pc, centre 
of the large number of individual and statistical distances 
reported in the literature) is 0.08 pc, i.e. 40 times shorter 
than the CII layer. An even more unfavorable situation 
occurs for a nebular distance of 980 pc, as adopted by 
Bianchi and Grewing (1987). 

We wish to stress here the great opportunities opened 
by the analysis of the zero-velocity columns in expand- 
ing nebulae (planetary nebulae, shells around Population 
I W-R stars, supernova remnants etc.): for the first time 
we can insulate a definite slide of nebula, thus removing 
the tiresome projective problems related to these extended 
objects. 

To be noticed the strict analogy between the zero- 
velocity column obtained from slit spectroscopy and the 
rest-frame image derived from Fabry-Perot interferome- 
try; in a certain sense, the rest-frame image is a zero- 
velocity column extended at all position angles (and this 
represents a considerable observational advantage). 



Moreover, the expansion velocity of a few objects (for 
instance the Crab nebula) is so large that the "zero- 
velocity column" analysis can be applied to very low res- 
olution data (slit spectroscopy or Fabry-Perot interferom- 
etry). 

Concerning our tomographic results, we confess a gen- 
uine satisfaction for the quality of these "preliminary" 
maps (clearly, the final goal of each tomographic study 
is a complete, three-dimensional model, but this will need 
a more detailed spectroscopic coverage) . 

Our figures confirm that the main body of NGC 40 
has an inhomogeneous, elongated barrel-shaped struc- 
ture ("an opened-up ellipse" following Mcllcma, 1995; "a 
slightly tilted cylindrical sleeve" following Balick et al., 
1987), with thin arcs emerging at both ends of the major 
axis. 

This morphology is normally interpreted as the result 
of the interaction of the fast stellar wind from the high 
temperature central star with the slow, inhomogeneous 
(i.e. denser in the equatorial plane) super-wind from the 
AGB progenitor (sec Balick et al., 1987, Mcllcma, 1995 
and references therein). 

Following Balick et al. (1987) the stellar wind cre- 
ates in NGC 40 an interior bubble at high temperature 
and ionization. The hot gas cools by expansion or radia- 
tive losses in forbidden lines when it interacts with the 
dense, external shell; [OIII] emissions arise in a relatively 
smooth, low density interface between the bubble gas and 
the Ha+[NII] filaments. 

As noticed by Meaburn et al. (1996), in this case the 
[OIII] image is expected to follow the Ha+[NII] one, and 
both are dominated by the swept-up shell. The alterna- 
tive explanation proposed by Meaburn et al. is based on 
Mcllcma's (1995) models: the nebular material is swept-up 
both by the fast wind and by the ionization front. So, in 
the first evolutionary phases, when the central star tem- 
perature is still relatively low, the two shells appear as 



10 



F. Sabbadin et al.: Tomography of the low excitation planetary nebula NGC 40 



distinct features: Ha and [Nil] prevail in the shell swept- 
up by the H ionization front, and [OIII] in the wind-swept 
shell. In the specific case of NGC 40, the second, faster 
shell has almost undertaken the first one. 

Which model is right? The answer provided by our 
data is ambiguous: on one side we obtain that the [OIII] 
and the Ha+[NII] emitting regions are disconnected (in 
many cases they are opposite, i.e. strong [OIII] is as- 
sociated with faint Ha+[NII], and vice versa); on the 
other side the density distribution across the shell never 
shows evidence of a double peaked structure implicit in 
Mellema's models. 

So, for the present, we can only say that: 

- in the equatorial region of NGC 40 the stellar wind 
is blocked by the dense nebular gas and [OIII] emissions 
occur in the innermost, faster expanding part of the "bell" 
density profile; 

- at the poles, where densities are lower than at the 
equator, the braking effect by the nebular gas is less effi- 
cient and [OIII] appears in extended zones, internally to 
the polar caps. 

A question arises: since we observe modest expansion 
velocities in the nebular gas, where (and how) is decel- 
erated the stellar wind, whose original speed is of about 
2000 km s- 1 ? 

It is evident that further, more accurate observations 
are needed not only for the nebula, but also for the central 
star. At the moment, mass-loss rates and terminal wind 
velocities reported in the literature for the WC8 nucleus of 
NGC 40 span in the range 3xl0~ 8 M Q yr" 1 and 2600 km 
s- 1 (Ccrruti-Sola and Pcrinotto, 1985) to lxl0~ 6 M Q yr" 1 
and 1800 km s _1 (Bianchi, 1992). Moreover, the presence 
of wind fluctuations was recently discovered by Balick et 
al. (1996) and Acker et al. (1997). 

The ionized nebular mass obtained from the observed 
H/3 flux (Sabbadin et al., 1987) results to be M ion = 2 - 
5 xlO -2 M (depending on the adopted values of F(H/3), 
c(H/3), distance, angular radius, electron density and elec- 
tron temperature), and a conservative upper limit to the 
nebular mass pumped up by the fast stellar wind can be 
put at lxl0~ 2 Mq. 

Even in the most favorable case (i. e. momentum con- 
servation), Bianchi's wind (10 -6 M Q yr -1 and 1800 km 
s _1 ) needs less than 20 years to produce the observed 
nebular acceleration (and in this case, where, and how, 
is dissipated the wind luminosity, amounting to 500 Lq , a 
value which is at least ten times the nebular luminosity?). 

A longer interaction time (500 years) is obtained with 
the mass-loss rate and wind velocity reported by Cerruti- 
Sola and Perinotto (1985); in this case the wind luminosity 
is 33 L Q , of the same order of magnitude of the nebular 
luminosity. 

A final question concerns NGC 40 as prototype of 
a small group of very low excitation planetary nebulae 
powered by "late" WR stars: are the peculiarities found 
in NGC 40 (in particular, the positive inwards velocity 



gradient and the radial chemical composition gradient) a 
common characteristic of these nebulae (e. g. SwSt 1, PM 
1-188, Cn 3-1, BD+30°3639, He 2-459, IRAS 21282, M 
4-18, He 2-99, He 2-113, He 2-142, He 3-1333, Pc 1-7 and 
K 2-16)? 

Though a definitive answer will need accurate imagery 
and spectroscopy of a suitable sample of objects, a first 
indication comes from the very recent results by Bryce 
and Mellema (1999), showing that BD+30°3639 expands 
faster in [OIII] than in [Nil]. 

9. Conclusions 

Long-slit echellograms at different position angles of the 
bright core of NGC 40 allowed to analyse the expansion 
velocity field (the higher is the ionization, the faster is the 
motion), the radial density distribution (a "bell" profile 
reaching Ae=4000 cm -3 ), the radial ionization structure 
(peculiar in the innermost parts), and the chemical abun- 
dances (a gradient is present, due to contamination of the 
hydrogen-depleted stellar wind). 

Moreover, tomographic maps are obtained, giving the 
spatial matter distribution along the cross-sections of neb- 
ula covered by the slit. 

All these observational results, discussed within the 
interacting winds model, point out our incomplete, "neb- 
ulous" knowledge of the physical processes shaping this 
peculiar planetary nebula where "the tail wags the dog" 
(as wittily suggested by Prof. Lawrence Aller, private com- 
munication) . 
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